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Abstract 



The functional annotation of genomes, construction of molecular networks and novel drug target identification, are important challenges that need 
to be addressed as a matter of great urgency 1 " 4 . Multiple complementary 'omics' approaches have provided clues as to the genetic risk factors 
and pathogenic mechanisms underlying numerous neurodegenerative diseases, but most findings still require functional validation 5 . For example, 
a recent genome wide association study for Parkinson's Disease (PD), identified many new loci as risk factors for the disease, but the underlying 
causative variant(s) or pathogenic mechanism is not known 6, 7 . As each associated region can contain several genes, the functional evaluation of 
each of the genes on phenotypes associated with the disease, using traditional cell biology techniques would take too long. 

There is also a need to understand the molecular networks that link genetic mutations to the phenotypes they cause. It is expected that disease 
phenotypes are the result of multiple interactions that have been disrupted. Reconstruction of these networks using traditional molecular methods 
would be time consuming. Moreover, network predictions from independent studies of individual components, the reductionism approach, will 
probably underestimate the network complexity 8 . This underestimation could, in part, explain the low success rate of drug approval due to 
undesirable or toxic side effects. Gaining a network perspective of disease related pathways using HT/HC cellular screening approaches, and 
identifying key nodes within these pathways, could lead to the identification of targets that are more suited for therapeutic intervention. 

High-throughput screening (HTS) is an ideal methodology to address these issues 9 " 12 , but traditional methods were one dimensional whole-well 
cell assays, that used simplistic readouts for complex biological processes. They were unable to simultaneously quantify the many phenotypes 
observed in neurodegenerative diseases such as axonal transport deficits or alterations in morphology properties 13, 14 . This approach could not 
be used to investigate the dynamic nature of cellular processes or pathogenic events that occur in a subset of cells. To quantify such features one 
has to move to multi-dimensional phenotypes termed high-content screening (HCS) 4, 15 " 17 . HCS is the cell-based quantification of several 
processes simultaneously, which provides a more detailed representation of the cellular response to various perturbations compared to HTS. 

HCS has many advantages over HTS 18, 19 , but conducting a high-throughput (HT)-high-content (HC) screen in neuronal models is problematic 
due to high cost, environmental variation and human error. In order to detect cellular responses on a 'phenomics' scale using HC imaging one has 
to reduce variation and error, while increasing sensitivity and reproducibility. 

Herein we describe a method to accurately and reliably conduct shRNA screens using automated cell culturing 20 and HC imaging in neuronal 
cellular models. We describe how we have used this methodology to identify modulators for one particular protein, DJ1 , which when mutated 
causes autosomal recessive parkinsonism 21 . 

Combining the versatility of HC imaging with HT methods, it is possible to accurately quantify a plethora of phenotypes. This could subsequently 
be utilized to advance our understanding of the genome, the pathways involved in disease pathogenesis as well as identify potential therapeutic 
targets. 



Video Link 



The video component of this article can be found at http://www.jove.com/video/3452/ 



Protocol 



1. Automated adherent cell culture process (Figure 1) 20 

1 . Prepare the automated cell culture system for the input of cell culture plates. Load consumables (e.g. pipette tips, cell culture plates, assay 
plates) into the system using the graphical user interface (GUI). Ensure there is sufficient cell culture media, phosphate buffered saline (PBS) 
and trypsin in the robotic system. 

2. Manually seed two omnitray plates with 2 x 10 6 cells per plate, of the SH-SY5Y neuroblastoma cell line. Maintain cells in Opti-MEM with 10% 
fetal bovine serum (FBS). Put the plates into the cell culture robot using the GUI. The cells will be incubated at 37°C and 5% CO2. 
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3. Choose which cell culture protocol needs to be initiated 20 . One can chose from an adherent cell culture process 22 , culturing and expansion of 
embryonic stem (ES) cells on mouse feeder cells 23 or culturing of suspension cells 22 . 

4. Select the adherent cell culture protocol (Figure 1) and ensure adherent cell line specific parameter files are adjusted so that the confluence 
threshold (the area of the omnitray plate that contains cells) is set at 70%. Set total trypsinization time to two minutes. 

5. Instruct the robot to prepare new omnitray plates, with a seeding cell number of 2 x 10 6 cells per plate. 

6. Input omnitray plates into the cell culture system using the automated adherent cell culture protocol (Figure 1). This protocol involves the 
following steps: plates are incubated and imaged until they reach the pre-defined confluence threshold. If cells do not reach the confluence 
threshold within 5 readings, plates are removed from the system. Upon reaching the user defined confluence threshold, cells are washed, 
trypsinized and counted. A pre-defined number of cells are added to new omnitray plates and if there are sufficient numbers of cells, a 
specified number of assay plates are transported to the deck and a defined number of cells dispensed into each well. Assay plates can be 
directly imaged using the integrated microscope or output from the system for further processing. 

7. Instruct the system to prepare 4 assay plates per omnitray plate, with a total of 5,000 cells per well. 

2. shRNA virus production and plating into assay plates (Time required: 6 days) 

1 . Grow bacterial glycerol stocks containing the shRNA vectors (Open Biosystems, TRC1) overnight in 2ml of Luria-Bertani medium media 
containing 100 ug/ml of ampicilin (Sigma-Aldrich). 

2. Extract plasmids following the manufacturer's protocol (Promega Wizard MagneSil Tfx). 

3. Produce virus using the RNAi Consortium High-Throughput Lentiviral Production (96 well plate) protocol 24 . Working with lentivirus is 
relatively safe because the virus particles used for transduction are replication-deficient and split-gene packaging strategies are used for their 
production. However, when working with lentivirus, additional biosafety procedures are necessary to minimize the risk to oneself and others 25 
All experiments must be conducted in a MLII or BSL2 safety level laboratory. All plastics (pipettes, plastic dishes, media) that have been in 
contact with lentivirus particles should be incubated with bleach for 24 hours prior to disposal. 

4. Calculate the multiplicity of infection (MOI) of the lentivirus by determining the percentage of GFP positive cells using the pLKO.1 GFP 
plasmid (Sigma-Aldrich). 

5. Plate lentivirus into assay plates, with an MOI of 3. 

3. Lentiviral transduction and neuronal differentiation of SH-SY5Y cells (Time required: 6 
days) 

1 . Cells are added to the assay plates (see step 1 .7). Load assay plates containing shRNA lentivirus into the automated cell culture system. 

2. After 24 hours, the media on the assay plates will be changed to Opti-MEM containing 0.5% FBS and 0.1 uM retinoic acid to begin the 
differentiation process. Differentiation of SH-SY5Y cells allows visualization of neuritic structures and synchronizes cell division. 

3. Continue incubation of the assay plates in the differentiation media for 5 days. This ensures maximal knockdown of target gene expression. 

4. On day 5, add 50 uM H2O2 to half of the assay plates for 24 hours to stimulate translocation of DJ1 to the mitochondria. 

5. On day 6, add Mitotracker CmxROS (Invitrogen) to the cells, at a final concentration of 200 nM per well and incubate at 37°C for 30 minutes. 

6. One can instruct the system to image plates directly using the HC imager or plates can be exported from the system for further processing. 

4. Automated immunostaining of assay plates (Time required: 2 days) 

Image quality is paramount for conducting a sensitive and reliable HCS. Damage to the cellular monolayer due to inaccurate pipetting can lead to 
poor image quality and irreproducible results. In order to minimize cell layer damage, the immunostaining was conducted using a robotic station. 
The procedure is similar to one that has been previously described 26 but has been customized to increase throughput and reduce consumable 
usage. 

1 . Fix cells with 100 ul of 4% paraformaldehyde pre-warmed to 37°C. Incubate for 20 minutes at room temperature. 

2. Wash cells with 200 ul of PBS for 5 minutes, 3 times. 

3. Incubate assay plates with 200 ul of PBS containing 0.1 % Triton (PBST) for 10 minutes. 

4. Wash cells with 200 ul of PBS for 5 minutes, 3 times. 

5. Incubate assay plates with 200 ul of block buffer (PBST with 5% FBS) for 1 hour at room temperature. 

6. Wash cells with 200 ul of PBS for 5 minutes, 3 times. 

7. Incubate with the following primary antibodies overnight at 4°C: 

1 . Goat DJ1 N20 (Santa Cruz, 5 ug/ml) 

2. Rabbit (3-III tubulin (Sigma-Aldrich, 1 ug/ml) 

8. On the following day, wash cells with 200 ul of PBS for 5 minutes, 3 times. 

9. Incubate assay plates with the following secondary antibodies 1 hour at room temperature: 

1 . AlexaFluor 488 donkey anti-goat (Invitrogen, 2 ug/ml) 

2. AlexaFluor 647 goat anti-rabbit (Invitrogen, 2 ug/ml 

10. Wash cells with 200 ul of PBS for 5 minutes, 3 times. 

11. Incubate cells with Hoechst (Invitrogen; 1 ug/ml) for 10 minutes. 

12. Wash cells with 200 ul PBS for 5 minutes, 3 times. 

13. Store plates at 4°C until they can be imaged. 

5. High content image acquisition and image analysis (Time required: 5 days) 

1 . Image a total of 30 fields per well using the 20x objective lens. Visualize DJ1 with the FITC filter set, the mitochondria with the TRITC filter 
set, (3-III tubulin with the Cy5 filter set and the nuclei using the UV filter set (Figure 3). 

2. Analyze the images using the Compartmental Analysis Bioapplication (Cellomics, ThermoFisher) to determine the average intensity of the 
Mitotracker signal within the mitochondria. (Figure 4B, F). 

3. To determine the average overlap coefficient between DJ1 and the mitochondria, analyze the images using the Cellomics Colocalisation 
bioapplication (Cellomics, ThermoFisher). Define regions of interest (ROI) as follows: ROI A - nucleus (Figure 4A, E), ROI B - mitochondria 
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(Figure 4 B, F). Exclude ROI A from ROI B to ensure analysis of only the cytoplasm. Define the mitochondria as target region I and DJ1 as 
target region II (Figure 4C, G). 

4. Analyze the images using the Neuronal Profiling bioapplication (Cellomics, Thermofisher) to trace the average lengths of the neurites from the 
(S-lll tubulin staining (Figure 4D, H). 

5. Image plates using the Opera LX automated confocal reader (Perkin-Elmer). Image a total of 30 fields per well using the 60x objective lens 
with water immersion. Visualize mitochondria with the 561 nM laser and nuclei with UV excitation. 

6. Analyze the images using the Spot-Edge-Ridge (SER) texture features algorithm. The SER-Ridge filter transmits intensity in pixels forming 
ridge-like patterns. The more fragmented the mitochondria, the higher the SER-Ridge score (Figure 8). 

6. Data normalization and analysis 

1 . Import the data from the image analysis software into the BioConductor CellHTS2 package for the R software environment (R version 2.11 .1 , 
BioConductor version 2.6). 

2. Logarithm base (2) transform the data prior to per plate median based normalization 27, 28 . Do not apply the variance adjustment per plate. 

3. To identify modifiers of a phenotype, use a two way ANOVA between the different treatment groups i.e. Scrambled infected untreated cells vs. 
scrambled toxin treated cells vs. target gene untreated cells vs. target gene treated cells (Figures 5-7). 



Mutations within DJ1 give rise to early onset-recessive parkinsonism 21 , but it is unclear how loss of DJ1 gives rise to the disease phenotype. It is 
known that cells deficient of DJ1 are more susceptible to oxidative stress-induced cell death and in response to oxidative stress, DJ1 translocates 
from the cytoplasm to the mitochondria 29, 30 . By constructing HC assays to monitor these phenotypes, we can identify genes that regulate or 
affect phenotypes associated with DJ1. This approach can help decipher the pathways within which DJ1 functions and that could be involved in 
disease pathogenesis. 

Example of an epistatic interaction with DJ1 (Figure 5): Knockdown of DJ1 in cells exposed to toxin results in a greater loss of cell viability 
(BAR-B: image-B) compared to cells infected with scrambled lentivirus (BAR-A: image-A). Knockdown of target gene A has a similar effect to that 
observed in cells with a DJ1 knockdown (BAR-C: image-C). Knockdown of both DJ1 and target gene A results in a significantly greater loss of cell 
viability than loss of either gene alone (BAR-D: image-D). This suggests an epistatic interaction between DJ1 and target gene A. 

Example of a gene regulating DJ1 translocation (Figure 6): When cells are exposed to a toxin, DJ1 translocates from the cytoplasm to the 
mitochondria, which is quantified by a higher overlap coefficient between DJ1 and the mitochondria (BAR-A: image A versus BAR-C: image C). In 
cells where target gene B has been silenced, less DJ1 translocates to the mitochondria when cells are exposed to the toxin. This suggests that 
target gene B is involved in the transport of DJ1 to the mitochondria. (BAR-B: image B and BAR-D: image D) 

Example of a gene involved in neuronal outgrowth (Figure 7): Knockdown of target gene C in wild type SH-SY5Y cells results in a significant 
increase in neurite length (BAR-B: image-B) compared to cells infected with lentivirus expressing scrambled shRNA (BAR-A: image A). This effect 
is lost in cells incubated with toxin (BAR-C and D). 

Example of a gene involved in mitochondrial morphology (Figure 8): Infection of wild type SH-SY5Y cells with shRNA targeting gene D results in 
a decrease in the mitochondrial SER-Ridge segmentation value (Figure 8, Image C and D) when compared to cells infected with scrambled 
lentivirus (Figure 8, Image A and B). 



7. 



Representative results 
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Flow chart for the cultivation of adherent cell lines 
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Figure 1. Outline of automated cell culture protocol. Please click here to see a large version of this figure. 
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Figure 2. Schematic overview of the screening process, image analysis and statistical methods used during the screening process: i) Cells are 
cultured until they are confluent and subsequently plated into the assay plates containing the shRNA lentivirus. Cells are differentiated for 5 days 
and toxin is then added to the plates for 24 hours. Assay plates are output from the system and immunostained. The amount of time taken for 
each of the processes is indicated in brackets, ii) Data is acquired using a HC imager (cell viability, protein translocation and neuronal outgrowth) 
and an automated confocal imager (mitochondrial morphology). Data are exported to CellHTS2 package within R, base (2) log transformed and 
normalized. Two-way ANOVA is used to identify significant interactions between the different variables. Please click here to see a large version of 
this figure. 




Figure 3. Composite images of cells acquired by HC imaging. A) Untreated cells, B) Cells treated with H2O2. DJ1 is labeled in green, the 
mitochondria in red and the nuclei in blue. Neurite staining is not highlighted. 
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Figure 4. Quantification of several cellular features obtained from a HCS. A-D) Untreated SH-SY5Y cells. E-H) H 2 0 2 treated SH-SY5Y cells. A, E) 
Nuclei segmentation and definition of ROI A; B, F) Identification and quantification of mitochondria, ROI B; C, G) Identification of DJ1, Target 
channel II; D, G) Identification and calculation of average neurite length. Cells close to the edge of the image are excluded from the analysis. 
Inset images are images prior to analysis. 
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Figure 5. Identification of a gene in epistasis with DJ1 (letters on the bars correspond to the lettering on the images). Images A through D are 
SH-SY5Y cells labeled with Mitotracker CMXRos (red) that was used for quantification of cell health. 
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Figure 6. Identification of a gene regulating DJ1 translocation (letters on the bars correspond to the lettering on the images). Images A through D 
are SH-SY5Y cells labeled for DJ1 (green), mitochondria (red) and nuclei (blue) that were used for quantification of DJ1 translocation to the 
mitochondria. 
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Figure 7. Identification of a gene involved in neuronal outgrowth (letters on the bars correspond to the lettering on the images). Images A through 
D are SH-SY5Y cells labeled for f3— III tubulin (green) and nuclei (blue) which were used for quantification of neurite length. 
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Discussion 



With the decreasing costs of HT/HC cellular screening systems, combined with the availability of powerful genome-wide tools to modify gene 
function, HT/HC screens are becoming commonplace in academia. The approach has already been successfully applied to diverse areas of 
research such as the identification of drug targets in cancer 9 - 31-33 and embryonic development 34-36 and even has potential for application in 
deciphering the pathways involved in neuropsychiatric disorders 37,38 . However implementation of such a system requires a significant investment 
of time and effort with process optimization often taking a minimum of 6 months. All steps, such as trypsinization times, pipetting speeds and 
seeding densities have to be adjusted, to ensure that cells are healthy and grow consistently. Prevention of bacterial contamination is one of the 
most difficult challenges facing automated cell culture with weekly cleaning protocols in combination with constant rinsing of all liquid bearing lines 
with 70% ethanol being necessary for contamination free cultures. It will be also necessary to improve the robotics so that additional instruments, 
such as confocal microscopes for higher resolution and -80°C freezers for compound storage can be integrated. 

There are also limitations that need to be addressed to improve the sensitive, speed and utility of this method to study gene networks and identify 
genes involved in pathogenic molecular pathways. 

To conduct a HT/HC screen and ensure that reliable data is collected, several aspects have to be optimized. First, reliability of the measurement 
is paramount and is dependent upon the robustness and sensitivity of the assay. For example, the assays described above are suitable for 
smaller screens, but are difficult to implement on a genome wide scale, due to the number processing steps required before image acquisition. 
Thus, one would have to construct stable cell lines expressing the reporter gene, which would allow for direct imaging and lead to decreased 
variation due to the reduced number of processing steps. At present, designing an assay that accurately depicts and reliably quantifies a 
phenotype of interest is a major bottleneck in the HC screening process. 

Many screens are conducted in mammalian cells using different RNAi libraries, all of which suffer from off-target effects, limited gene silencing 
efficiency and incomplete genome coverage. Thus libraries need to be made which are more specific, potent and have better coverage. Efforts 
are underway to create such libraries (http://www.sigmaaldrich.com/life-science/functional-genomics-and-rnai/shrna/library-information.html) and 
it is hoped such endeavors will improve the reproducibility of HT/HC screen hits. 

A limitation of many large scale cell based screens are that they are conducted in neuroblastoma cells because they can be genetically 
manipulated and cultured to large numbers with relative ease. However, the relevance of 'hits' identified in ex vivo cell culture models to in vivo 
function is questionable, especially as the brain consists of highly specialized cell types that form a dense and complicated network of synaptic 
connections to function as a highly integrated unit. As a consequence, it is common that hits identified using the screening approach described 
above, are validated in secondary screens using additional techniques and in more physiological relevant models 39 . To improve the translation of 
hits identified during HCS, more representative and sophisticated models, such as primary cells and differentiated stem cells or co-culture 
systems need to be developed and adapted for HT/HC approaches. 

With a combination of automated cell culturing and HC imaging one can rapidly gain new insights into how neurons function and determine which 
pathways are important to disease development. Combining HCS/HTS data with information generated from other 'omics' approaches, it will then 
be possible to construct a systems biology overview of brain diseases, thus facilitating therapeutic development. 
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